Overlap of QRPA states based on ground states of different nuclei 
—mathematical properties and test calculations— 
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The overlap of the excited states in quasiparticle random-phase approximation (QRPA) is calcu- 
lated in order to simulate the overlap of the intermediate nuclear states of the double-beta decay. 
Our basic idea is to use the like-particle QRPA with the aid of the closure approximation and calcu- 
late the overlap as rigorously as possible by making use of the explicit equation of the QRPA ground 
state. The formulation is shown in detail, and the mathematical properties of the overlap matrix 
are investigated. Two test calculations are performed for relatively light nuclei with the Skyrme 
and volume delta-pairing energy functionals. The validity of the truncations used in the calculation 
is examined and confirmed. 

PACS numbers: 21.60.Jz, 23.40.Hc 



I. INTRODUCTION 



One of features of neutrino physics is its interdisci- 
plinarity. Neutrino physics is important for a better un- 
derstanding of particle physics in terms of lepton- number 
violation, Majorana nature of neutrino, and neutrino 
mass, e.g. Refs. [lrl!, (there are many textbooks eluci- 
dating neutrino; see for example Ref. Q), which are as- 
pects of particle physics beyond the scope of the standard 
model. Neutrino physics is also very interesting from the 
viewpoint of nuclear physics [5l-fl3| . One of the few meth- 
ods used to determine the neutrino mass requires accu- 
rate calculations of the nuclear matrix elements in the 
neutrino-less double-beta (OvBB) decay along with its ex- 
perimental half life, e.g. Refs. For details on other 



methods to determine the neutrino mass, see e.g. Ref. [14 1 
(shape of the /3 decay spectra and other particle-physical 
methods) and Ref. [15( (cosmological method). The pri- 
mary task is to determine the neutrino mass accurately, 
on the other hand, this is a very good and challenging 
opportunity for theoretical nuclear physics to test if the 
techniques developed so far in this field are useful for 
solving the problem of other field. This is particularly 
because most of the nuclei providing the ground of the 
Qvfifi decay are heavy nuclei, and therefore many-body 
correlations have to be taken into account along with 
large wave-function space. 

In this study, we take the first step 1 in the attempt to 
calculate the nuclear matrix elements of the Qv/3/3 decay 
by making use of a method different from the traditional 
ones; First, we use the like-particle quasiparticle random- 
phase approximation (QRPA) pjj (formulation of the 
QRPA) [lj| (application to the CV/3/3 decay suggested) 
for axially-deformed nuclei, which can be applied after 
the closure approximation is used. This approximation 



1 Parts of the formulation and the test calculations in this paper 
are reported in Ref. |l6t . 



has been proven to be good in the 0^/3/3 decay by the 
analytical argument [ll| and several realistic calculations 

HE [312. 

In the application of the QRPA to the 0^/3/3 decay, two 
QRPA-state spaces are obtained via calculations based 
on the initial and final states of the decay, and the prod- 
uct of the two projection operators to the QRPA-state 
spaces is inserted to the middle of the two-body Of/3/3 
transition operator. Secondly, in our approach, the over- 
lap of the intermediate states obtained by the two QRPA 
calculations is calculated more accurately than ever. A 
simple approximation and a few variants [H, [23l - [26| have 
been used for calculating the overlap. The importance of 
the overlap of the intermediate states is pointed out in 
Ref. [23[ in terms of deformation. It is reasonable that 
the overlap is sensitive to the difference in the deforma- 
tion of the initial and final states; this raises a question if 
the differences in other properties affect the overlap. We 
can address this question comprehensively by treating 
the ground-state wave function of the QRPA explicitly, 
and here we demonstrate the feasibility of that treatment 
and investigate mathematical properties of the overlap. 
The equation of the QRPA ground state has been known 
for decades, e.g. Ref. (27j; however, to the best of our 
knowledge, our study is the first that carries out the in- 
volved numerical calculation rigorously. Few researchers 
have attempted to calculate explicitly the QRPA ground 
state in subjects other than the study of the Of/3/3 de- 
cay [27M29I] , and they have mostly used crude approxima- 
tions. The probable reason for the rarity in attempting 
this calculation is that it is possible to obtain the transi- 
tion strength of the QRPA without treating the explicit 
ground-state wave function (l7j . In this light, neutrino 
physics provides further motivation to develop techniques 
of nuclear theory. 

The third feature of our approach is the use of the 
Skyrme energy density functional (30T [3l| . It is of in- 
terest from the viewpoint of nuclear theory to investi- 
gate how a phenomenological approach developed so as 
to reproduce as many experimental data as possible in- 
cluding the masses and the root-mean-square radii of the 
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ground states can be successful in describing other nu- 
clear properties. The Skyrme energy density functional 
has been used for providing the Hartree-Fock field to 
the calculation of the nuclear matrix elements [HJ . We 
use the Skyrme-plus-pairing energy density functional to 
solve the Hartree-Fock-Bogoliubov (HFB) [TtJ and the 
subsequent QRPA equations self-consistently. The self- 
consistency assists in strengthening the reliability of the 
calculation. 

The standard method in the category of the QRPA for 
calculating the nuclear matrix elements is the proton- 
neutron (pn) QRPA [33| . It has been argued that the 
Pauli correction terms are necessary to include in the 
calculation of the intermediate states. The renormal- 
ized pn-QRPA [33 - [3t| has been used for including the 
Pauli correction terms, and later the Ikeda sum rule [37j 
was satisfied upon using the fully renormalizcd pn-QRPA 
[Uliil. The self-consistent HFB and pn-QRPA calcula- 
tions have been performed in Ref. poj . The importance 
of the particle-particle interaction has been pointed out 
in Ref. j4lj , and subsequently the proton- neutron pairing 
correlations have also been included in Ref. [42[. Fur- 
ther, a pn-QRPA calculation using the unitary corre- 
lation operator method has been performed for taking 
into account the short-range correlations [43[. As previ- 
ously mentioned, another improvement as regards the pn- 
QRPA is its extension to deformed states Q EE)]. Thus, 
the pn-QRPA has been improved up to a very advanced 
level in the past few decades. Nevertheless, as is well 
known, the problem of the systematic difference in the 
nuclear matrix elements between different approaches has 
not thus far been resolved [1, EH; in particular, there is 
a difference of a factor of two between the pn-QRPA and 
the shell-model approach. As for approaches other than 
the pn-QRPA, these can be found in e.g, Refs. [47l - [5l| 
(the shell-model), [52| (t he projected HFB), [|53[ (the mi- 
croscopic IBM), and [54[ (the energy-density functional- 
plus- generator-coordinate method) . 

One of the advantages of our method is that the feasi- 
bility of the like-particle QRPA calculation is fairly high 
for any nuclei except for the transitional ones between 
the spherical and deformed regions. To the best of our 
knowledge, the collapse of the like-particle QRPA due to 
the pairing fluctuation does not occur [55| as long as the 
strength of the pairing energy functional is determined 
so as to reproduce the pairing gaps obtained from ex- 
perimental odd-even mass differences (56|. Another ad- 
vantage is that the calculation of even-even nuclei is free 
from a problem that the last odd particle may not be 
approximated very well by the HFB calculation, if the 
coupling of the last particle to the nucleus is weak. The 
drawback of our approach is that since the closure ap- 
proximation is not good for the 2vj3j3 decay, e.g. [13j |. 
the reliability of our method is difficult to prove by itself. 
Perhaps it is necessary to rely on other methods to obtain 
a reference value of the representative energy of the in- 
termediate states which makes the closure approximation 
exact. If the effects of the higher-order many-body cor- 



relations beyond the QRPA arc minor, and a sufficiently 
large space of the intermediate states is used, then the 
question is whether the pn-QRPA and the like-particle 
QRPA provide similar nuclear matrix elements. The an- 
swer is not trivial, because the many-body correlations 
treated in the two QRPA methods are different. Thus, it 
is worthy to compare the two methods numerically. 

This paper is organized as follows: Section HT1 presents 
our basic scheme to calculate the nuclear matrix elements 
and the detailed formulation for calculating the overlap 
matrix elements of the intermediate states. The ana- 
lytical properties of the overlap matrix are discussed in 
Sec. Mil for simplified cases. Section ITVl provides techni- 
cal information regarding the numerical calculations for 
this paper. The calculated overlap matrix elements are 
shown in Sec. |Vl and the truncation approximations are 
examined in detail. Section I VII summarizes the study. 



II. FORMULATION 

The axial and parity symmetries of the nuclei are as- 
sumed throughout this paper. The z-component of the 
angular momentum is denoted by j* for nucleon state 
a and by K m for nuclear state m. The terms 7r Q and 
7r m are used to indicate the parity. All of the matrix el- 
ements used in the numerical calculations of this paper 
are real, although it is not assumed in the formulation of 
this section. Hereafter, we call the like-particle QRPA as 
simply the QRPA. For complete equations of the nuclear 
matrix elements, see e.g. Refs. UdElI. 



A. Application of QRPA to our method 

As mentioned in Sec. HI we assume that the closure 
approximation has been applied to the nuclear matrix el- 
ements of the 0is/3(3 decay. One of the components of the 
nuclear matrix element arises from the double Gamow- 
Teller operator and this component is now written 
as 

Mgt = (F\ E h + {nj,E a )a{i) ■ a(j)r + (i)r + (j)\I) 

ij 

a/3 a'f3' 

xt + {1)t + {2)\P'P){F\cIcI W \I), (1) 

where \F) and \I) denote the final and initial nuclear 
states of the decay, and the ground states of the QRPA 
are used. h + (rij,E a ) is the neutrino potential Q with 
fij = \fi — fj\, and E a being the average energy of 
the intermediate nuclear states. The index i (j) indi- 
cates a nucleon, <r denotes the spin-Pauli matrix, and 
t + denotes the raising operator of the z-component of 
the isospin. An arbitrary single-particle basis {a} is in- 
troduced, and the creation and annihilation operators of 
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the single-particle state are denoted by c a and c a , re- 
spectively. 

We introduce the creation and annihilation operators 
O 1 ^ and 0^, respectively, of the excited state m of the 
QRPA based on the initial state, and those based on the 
final state are denoted by O^} and 0^ n . The same kind 
of index m is used for specifying the QRPA states based 
on the initial or the final state. The states \F) and \I) 
are defined in the QRPA by 



oL\i) = o, 

O£\F)=0. 



(2) 
(3) 



Inserting the product of the two completeness equations 
of the relevant space 

i = \i)(i\ + J2oi!\i)(i\o! n 

m 

+ E oilo I „l\i){i\o I m2 o I mi + -- - , 

mi77l2 

1 = \F)(F\+J20^\F)(F\OF 

m 

+ E O^Of r l\F)(F\0^ 2 Of ni + ■ ■ ■ , (4) 



to the middle of the product of the single-particle oper- 
ators, we get 



(Ou) 
GT 



EEv^'IM^aMlWCS) 

a/3 a'P' 

xr + (l)r + (2)|/3/3')E^l c « c ^°™ t l^) 

mm' 

x(F\0^0 I Jl t \I)(I\0 I mt cp l ep\I). (5) 



We assume that the higher-order terms with respect to 
Om or Om do not have contribution to the two- 
particle transfer matrix clement, for example 



(I\O* ni O I m9 cp,cp\I)=0. 



(6) 



(In the QRPA order, this equation holds exactly.) The 
nuclear states 0^\F) and 0% 1-0 in Eq. © are the inter- 
mediate states mentioned, and the overlap of these two 
intermediate states is not equal to 5 mm i in the QRPA. 
The conditions for the product 

{F\ci^,0^\F){F\0i02 / \I){I\0 I m ,cp,c p \I) 1 (7) 

in Eq. ([5]) to be finite are 



3a + 3a' =3/3+ 3/3' 

for the single-particle states, and 



n m n m Ja Ja'l 



(8) 



(9) 



for the intermediate states. For an arbitrary pair of a 
and a', there exist (3 and /?' satisfying condition (JSj) and 
the condition that the two-body matrix element of the 
double Gamow- Teller operator is finite. Thus, the QRPA 
solutions are necessary for all (K m ir m ) for which Eq. ([5]) 
is convergent; in other words, there is no selection rule 
for the intermediate states. 



B. Formulation of overlap of intermediate states 

In this subsection, we show the detailed equations for 
calculating the overlap matrix elements 



F {m\m') I = {F\Ol l O I 7 l,\I). 

Hereafter, we use the simplified notations K = K Tl 
7T = 7r TO . We express \I) and \F) in the form (271 ], 



\F) 



^ n ex p 

K it 

i 



Vj 



I cxp 



K' 



l/>, 



(10) 

and 

(11) 
(12) 



where \i) and |/) denote the HFB ground states of the 

nuclei described by \I) and \F), respectively, and v\ K n ' 

and Vp 71 ' denote the generators of the QRPA ground 
states. The terms JVi and J\fp indicate the normaliza- 
tion factors. We have [0^,0 TO '] = in the QRPA if 

(Ktt) =^ (K m iTr m i), and hence, Vj W ''s and t)^ "''s 
with different values of (K'tt') arc determined separately 
by using 



O^, exp 
C exp 



-(•K'm"W) 



\i)- 
\f) 



0, 



0. 



(13) 
(14) 



General quasiparticle bases, which are not necessarily 
the diagonal representation of the HFB Hamiltonian, are 
introduced by using \i) and |/} as the vacuum state, that 
is, 



<l/) = o, 



(15) 
(16) 



where // = (q^, 7r M , j* i M ) denotes the label of a general 
quasiparticle state. The term indicates proton or neu- 
tron, and ifj_ denotes a label specifying the general quasi- 
particle state in the subspace (g^, 7r M , j*). The notation 
— \i is used for expressing (<j M , 7r M , — j* i^). The genera- 
tors Vj K 71 •* and Vp 71 ^ are written as 



v 



(K'tt') 



E (17) 



,(ifV) _ 



^HV^'vi a n a u V V • 
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It is to be noted that a T Ja^ and a^aj, in Eq. (fl~7|) arc 
the fermion images of bosons. In relation to this, we in- 
troduce the condition that Cj^ Z„? docs not vanish only 

ftjp+Hl = K ', +Ju' = - K '> and 7r M 7r ^ = Kp-'Ku' = 7r'. 

We order the general quasiparticle states and place the 
restrictions of pi < v, p! < v' in Cj^Z^f without losing 

generality. These conditions are also applied to Zl* ' • 
If K' is equal to 0, Eqs. (fTTj) and (fT5|) contain the same 
product of the creation operators twice; that is our choice 
of convention for simplicity of the programming of the 
code. 

Solving the QRPA equation, we obtain 



o Jt , 

m 

o F ] 

77V 

o F , 



(\r Irri' 
pis 

( vim' * I I 



Y Im a 1 a 1 



E 

p<v 

E 
E 



El v Fm'* F F xrFm'* Ff Ff 
[X nv a v a^ -Y_ fl _ v a_' fi a_l 



yFm'FiFi v Fm' n F n F 



(19) 

(20) 
,(21) 
, (22) 



where j* + j z v = K m < and -k u -k v = n ri 
We define matrices 



x (K'-k')I 

y {k'-k')i . 





r (K' v ')I 
°ll,-n-n' 


r (K'ir')I 


c (K'ir')I 
nn' ,—ri—n 


Y n 
JL n 




yll 


x%) ' 


1 -i-i 


• Y IM \ 

■ ) • 


Y n , ■■ 

— n—n 


• Y IM , / 

—n—n / 



,(23) 

(24) 
(25) 



where the QRPA solutions of the (K'ir') are used. Matri- 
ces C( K '^ F , JSf(* rV ) F , and y(^ ,7r ')F are a i so introduced 
in the same manner. 

Subsequently, we obtain C^ K w ' J and C^ K 71 ^ F ignor- 
ing the exchange terms (the quasi-boson approximation 



Here, I denotes the unit matrix with the size of matrix 
D, and D nv is not equal to zero only for those \x and 
v satisfying j* + j% = and ir n ir u — +. We place the 
restriction of j* > in Eq. (f2"T)) . The unitary transfor- 
mation from the basis {a F a^A to the basis {a^",a^} 
is given by 



= Efe 1 v t +«< 



(30) 



Y,{nr<-+n F ->-i) . (3D 



with = j*, and tt h = 7r M ' . The matrix elements of the 
unitary transformation can be calculated as 



rpIFl 



T, 



IF2 



d 3 rJ2(U^(r,a)Ul(r,a) 
+V F *(r,a)Vl(r,a)) , 

+U F n ,(r,a)V^r,a)) , 



(32) 



(33) 



by usingthe wave functions of the general quasiparticle, 
see e.g. [581 ] . 



V.Ur,-a) 



(34) 



and those associated with the state F, where a = ±1/2 
is the z-component of the spin. Z7 M -„ is obtained from 



D 



( _}_rpiF2 

I J<IFl 



(35) 



q{K'-k')F _ 



1 



1 + 5kq 
1 



y-(*rV)J 



1 



1 



?A'0 



r (K'7r')i 



1 



X(K'n')F 



(26) 



where the suffix T indicates the transpose of matrix, and 
it is assumed that 1/X^ V ) J and l/X ( . K ' n '^ F do not 
have a singularity. 

The relation between the two HFB states can be writ- 
ten as, see e.g. [ll 



|l) = M exp 




(27) 
(28) 

(29) 



where the matrices used are defined as 



T 



IFl 



T 



IFl 



rpIFl 
1 11 



rpIFl 

1 nl 

1 1-1 



rpIF2 

1 n-l 



IFl 



rpIFl 

nn 



rpIF2 

1 l-n 



T, 



IF2 



(36) 



(37) 



and it is assumed that 1 /T docs not have a singularity. 
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Now, we expand and truncate the overlap matrix element with respect to ip^ * ^ and Vj K 71 ^ 



(F\OiO%\I) = j^rjjrifl II ex P 

1 F KflTt 



.(Kl7Tl)t 



ozo% n ex P 



Ko.TT-2 



m \{f\o F m o%\i)+ J2 ((/l^ 17ri)+ o^o^K) + (/|o^o^4 Kl7ri) K) 



K"l7Tl 



Kl7Tl 

X = 



7v/7V> ' 



■A/j 



Ki7ri 



i 







cxchj 



~,(^l*-l)t-(-Kl*-l)l 



Vcxch 



, p(Ki7ri)J _ ^(KiTri)! p(Ki7ri)I _ p(Ki7ri)l\ 
{M^vf^f) = (/|4* 17ri)t 4* 17ri) |/)bo S on + 17ri3t 4 Xl7ri) |/>exch, 



(ifiTrOt-fifiTi-x) 



(/|fif^Ot^*0| /)(otdi 



/)bo SO „ = (1 + ^ l0 )Tr (c(tfi-i)^C'(tfi-i)^ ; 



/, , r iWMM^/ ^(ifiTrOF ^(_ffi7ri)F ^(KTittiJF ^-.(KiTTi; 



_ ( j{K 1 tt 1 )F _ ^(Jfi7ri)F £,(iCi7ri)F _ ^(ifiTnjF 



The fourth-order terms in Eqs. (|40|) and (|41|) are approximated by the following quasi-boson terms 



1 

4 



^ | (2+ 6<5 Kl o) (Ti^C^^^C^ 1 " 1 '^)) 2 + (2 + 14^ l0 )Tr(C7 (Kl7ri)/ C ( ^ l7ri)7t ) 2 | : 



Klin 



i/) 



J^{(2 + 6<fe l0 ) (Tri&x^Ctx^W)) 2 + (2 + U^Tr^^C^ 1 ^) 5 



(38) 

(39) 

(40) 

(41) 

(42) 
(43) 



(44) 

(45) 
(46) 



(47) 



(48) 



(49) 



We test up to the second-order terms 



with respect to v^ 1 ™ 1 ^ and 0j 27r2 ^ but only with 



(/I 



(50) 
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(A'itti) = (K 2 ir 2 ) in Eq. ([38j) 2 . The terms up to the 
fourth-order arc included in the normalization factors Mi 
and Mf, because its convergence is slow with respect to 
the {i-expansion compared to the un-normalized overlap 
matrix elements, that is, Eq. (|38|) without Ml. The rea- 
son for this difference is that the bra and ket HFB ground 
states are identical in the normalization factors, while 
these states are quite different around the Fermi sur- 
face in the un-normalized overlap. Due to this difference, 

I 

The first term of Eq. (|38| is obtained 



The se 



The third 



Further, we can write the fourth term of Eq. p8|) as follows: 



high-energy excitations leaving the configuration around 
the Fermi surface intact do not significantly contribute 
to the un-normalized overlap matrix. Hence, the un- 
normalized overlap has less of the major terms than the 
normalization factors, and Vp 1 ™ («j 1 ) in the un- 
normalized overlap has a smaller effect than in the nor- 
malization factors (see the numerical results in Sec. IV)) . 



(51) 



(52) 



(53) 



li<v n' <i>' 



id term of Eq. (f3"8")) reads 

Ml E (f\*?™»o*(%.\i) 



KxTTX 

= ^E E E E c ^u' F * x ^ x lTu 2 (/i^af <<<<M*> 

K\TT\ flUfi'v' //l<fl M2<^2 

ka yFm* v-Im' ( ^(.KiTijF* _ ~,(KxTrx)F* 

J / j / j / j / j — fix— Vl H2"2\ — vx— /J.x,ixv ^—nx — vx^v 

K\TT\ flV Hx<VX [12<V2 

r (Kxnx)F* _ r (Kx-Kx)F* r (K 1 Ti 1 )F* _ r {Km 1 )F* ^(KxirxJF* r (K 1 -K 1 )F* 
fj,u,— Ux — ^i ^v. — fii—vi ' — Vxv, — MiM — Mi^ — ^iM — ^l^.M — Mi ' —Mi v,il — v\ 



r (K 1 TT 1 )F* _ r (Kxwx)F* _ r (Kx-Kx)F* r (Kxwx)F* \/f\ n FFIiIU-\ 
"l.— M— \ix,—vxv M~ v\ ,v—\l\ ' "-"/i — Vl I \J I U M V M2 V2 I / 



term of Eq. ([38]) is given by 
M E (f\OiO%4 K ^\i) 

= M E E E EE^^tti/la^^a^^^^Jti,) 

K\-K\ fl<V /j'<I/' Ml "1 M2^2 

-WV VV \" x Fm* Y Im' f _ r (Kxnx)I r (K l -K l )I 

z_^i ^ v 1 -n'-v'\ ^nxfi 2 ,-i>'-fi' ' ^fixH2,-n'-v' 

K\-k\ i±<v \x' <v' M1M2 

^{Kx-kx)! ^(Kx-nx)! _ ^(Kx^x)! ^(Kx^x)! , MKxtvx)I _ ^(XiTn)/ 

— V 1 — fJ, f ,(J,l(J,2 — /X 7 — I/ y ,jUljU2 \l\—[l' % \l% — v' fll —f' .(J.2— fl' [l\—\x'-, — v'\l% /Xi— f',— jU'jLi2 

r (A'i7ri)J r {Kx-Kx)I r {Kxitx)I r {Kx-Kx)I \(f\n F n F n I ^n I Ui\ 



i<"l7Tl 



^mm' * / y / , / y ^^u,fi'v' ' tl3V3,H&Vi / , / , Ml^l M2^2 

KxTTX fLV\l' v' M3^3M4^4 Ml<^lM2<^2 



jn2 . f (^i7Ti)X2 ylm j _ p(Kx-Kx)I , ^(ifiTri)/ 

mm' y / ; / y / < HI ^1 U )X1 t^l / j 1 —H2—V2 \ M3^3,— ^2-^2 ' ftS— f-2,V3— "2 

-ftTi7Ti Ml^l M2<f2 

M3-M2,-^2f3 ' -M2^3,M3-f2 — M2 1^, — "2^3 M3-"2,f3-M2 M3 — "2 1 ~ /*2 "3 -^2^3:M3-M2 

^(KlTTi)/ _ ^(ifiTTi)/ \ /r fi \ 

— "2^3, — M2M3 -M2— 1^2,M3^3 i ' V UU 7 
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p3 . _«y Ylm'p(K 1 TT 1 )X3\r^YFmffi(Knr 1 )F*_ ( ^(K 1 TT 1 )F* 

mm' J / j / j / j F.ft^ / j I fl'v' fiv,fj.' v' \ip,' ,vv' 

K\7t\ 11V /i2<^2 ll'v' 

^{Km^F* _ c ,(K 1 ti 1 )F* ( ~,(K 1 -K 1 )F* p(Ki7Ti)F* _ ^{Knr^F „{Kiirx)F* 

' [Ifl' ,V f V li'l/illU' fl' fl ' III/' ,Ufl' til/' ,fl' V ' V 1 

3 = e E ^c&</tf«M^taWti<> > ( 59 ) 



L< = ^EEEEE( 



jUI/, — [i 1 — v 1 /i— ii',i/— i/' im—ii',—v'v 



Kim p'v' ftv ^31/3 

_ r (K 1 7t 1 )F* _ c ,(K 1 tt 1 )F* r (K 1 n 1 )F* _ r (K 1 ix 1 )F* _ MK-Lit x )F* , MK-Liri)F* MKm^F* 
—fi'H^v—v' — ju'/z,— v'v ' /*— i/',i/— /a' /*— i/',— /i'i/ — Z/'/i,!/— /a' ' — it'll,— p,'v ' — /a' — i/',/AI/ 



A \J \ U, v "li >3 "3 II I A»3^3|— /*4 — ^4 /J3-M4, 1^3-^4 " T /i3~ A*4; — ^4"3 ' — ^4 M3 ,^3 - ^4 M4 M3 , - ^4 K3 



, fi(Kiiri)I _ ^(Kx-nx)! _ ^(K^^I ^~i(Kitti)I ^(Kxitx)! \yFm ylm rnrs\ 

' A»3— Vi,V3— lii Vi,—IHV 3 — 1/4^3,1/3 — ^14 L ^— I/ 4( L4 3 ,— li 4 l/ 3 ' — M4 — "4 ,M3"3 J — it' — 1/' 1 — (U4 — ^4 ' V 0U -' 

The generalized expectation value of the multiple fermion operators can be calculated by using the generalized 
Wick's theorem In particular, that of four operators can be written explicitly 

,F r,F_J\ r,l\U\ - 1 ( lt\„F „Fu\lt\„I\„I\\ A \ _ lt\ n F. n l\\i\lt\ n F. n I 1 \\i\ -L I f U*! „'t U\ I f \ n F n I\ I 



(/K'<< 2 <|i) - 7^) ((/K'<'N)(/I« T H> - </K< \i)(fK'<l\i) + (fK<l\WK>a£\i)) , (6i) 

with the following contractions 

^ (62) 

-<J>^<H) j;<o, 

E T i-W^4 - E E T l F l rt ^-v (f\a*<£p> |*> , £ > , 
n^o-tli) = <! 4 (63) 

J_ t ifi _ J_ V T IF2 ,D , v- > 
! ^ - (64) 



Here, t M is a phase arising from the time reversal of a III. ANALYTICAL PROPERTIES OF 

fermion state as OVERLAP MATRIX 



A. Simple model 



T\a}= t* a \-a), (65) 



Let us discuss the simple model shown in Fig. [T] for 
investigating the analytical properties of the overlap ma- 
trix elements of the QRPA states. A single-particle basis 
where T denotes the time- reversal operator [5||. is shared by the initial and final states of the 0^/3/3 de- 
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A' 




\f) 



10 



FIG. 1: (Color online) Simple model of final and initial 
states of OvPfi decay. The doubly-degenerated levels are in 
the relation of the time reversal. Levels A and A' are referred 
to in text. The notation p (n) denotes protons (neutrons). 



cay with no pairing field. In this model, we assume that 
\F) = |/) and \I) = \i). When the QRPA is applied to 
this model, only three types of excitations are possible: 
two-particle addition, removal, or one-particle-one-hole 
excitation. Then, only two overlap matrix elements of 
the excited states are finite; one is the excited state with 
both levels A and A' (see Fig.[T} occupied, and another is 
the excited state with none of those levels occupied. All 
the other overlap matrix elements vanish. This feature 
can be quantified by the following measure 



5 = Tr(G t G)/dim(G), 



(66) 



where G denotes the overlap matrix, and dim(G) denotes 
the dimension of matrix G. If G is unitary, then, of 
course, S is equal to 1. The S value of our simple model 
is 2/dinuS, which is of order of 10~ 4 or smaller with the 
dimension of the realistic calculations. Therefore, the 
overlap matrix discussed is highly non-unitary. One of 
its implications is that O m cannot be represented by a 
linear combination of 0\ n , and O^, . This is also seen 
from the nature of the Bogoliubov transformation. Since 
aF} and are represented by a linear combination of 



J1 



and a u i , O m includes the bilinear term 



Jlni. 



The 



appearance of this term is certain, because the two nuclei 
considered have different configurations. And, ajja 1 ^, is 

bilinear with respect to Olj, and O m ,, according to the 
boson-expansion theories pj], [60[ . 

The above argument using S indicates that the overlap 
matrix of the QRPA is not close to the one obtained from 
the exact many-body states at all; the exact one has the 
absolute value of every diagonal matrix element equal to 
1. We need to recall the nuclear matrix element [Eq. ([5])] 
and Eq. ((6|) in order to understand the implication of this 
mathematical property of the overlap matrix. The inclu- 
sion of higher-order excited states such as O^J O^J ]/) is 
necessary for having many diagonal overlap matrix ele- 
ments close to 1; however, these states do not contribute 
to the two-particle transfer matrix elements. Therefore, 
the QRPA has the possibility of being an efficient ap- 
proximation to the nuclear matrix elements irrespective 
of the deviation of the overlap matrix from that of the 
exact many-body states. 



The non-unitarity of the overlap matrix is the reason 
why we do not use the boson representation in the cal- 
culation. In fact, we have developed a code to use the 
boson representation disregarding the non-linear terms 
of the transformation between the two boson bases and 
re-orthonormalizing the transformation. Consequently, 
the absolute values of some overlap matrix elements ex- 
ceeded one by more than an order of magnitude. Thus, 
this artificial unitarization method using the boson rep- 
resentation cannot be accepted. 



B. Identical initial and final states 



We assume that |/) = \i) in this subsection. In this 



case, v 



(iCV)t 



[Eq. (fTT)) ] is equal to v' expressed as 

-t -at -f>t 

V =2^ V rn" V rn"> 

m " 
fi'v' 

*t f - = E ( jc 



(JLV 



(67) 
(68) 

(69) 



where the suffixes / and F as well as (K'n') are omitted. 
The condition K' ^ is also assumed; however, this 
assumption is not essential. It can be shown in the QRPA 
order that 



m" ' m 



8m' 



m" in' • 



(70) 



Let us suppose that the backward amplitudes of the 
QRPA solutions m and m' are very small. Using Eq. ([70]) 
and ignoring (i\0\ ■ ■ ■ 0\ n \i) ~ Y n with n > 1. we obtain 

/<m|m')/ ~ ^(i\O m Ol,0- + ^)(l + v)\i) 

~ ^(i\O m Ol,\i)(i\{l + v^(l + v)\i) 

= 5 m m'- (71) 

It is assumed in deriving the last expression that the nor- 
malization factor J\f 2 is calculated up to the same order 
with respect to v as that of the denominator. This deriva- 
tion implies that the truncation of the w-expansion does 
not affect the overlap matrix elements of the QRPA solu- 
tions that do not have the backward amplitudes. Thus, 
the non-collective states are expected to satisfy Eq. (|7T|) 
fairly accurately. From our numerical calculations shown 
in Sec. |Vl we confirmed that this expectation was correct 
for all of the non-collective states with a deviation less 
than 10 -4 . The deviation from Eq. (fTT]) with m = m' 
of the relatively collective real states among the QRPA 
solutions is around 0.01, and that of the two spurious 
states associated with the particle number is around 0.5. 
From this deviation, our method should be applied only 
to the cases for which the break in the particle number 
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conservation is not so large that a large deviation from 
Eq. (fTTj) for the spurious states does not significantly af- 
fect the nuclear matrix elements that we finally require. 



IV. DETAILS OF CALCULATIONS 
A. Procedure 

We use the code of the HFB approximation explained 
in Refs. [6ll - [o1| and that of the QRPA developed by us 
fall ]. Although the two codes have been developed in- 
dependently, in both codes, the wave functions are ex- 
pressed in the cylindrical coordinates. The wave func- 
tions are interpolated by the B-splines, see e.g. [65l f66l] . 
and contained in a cylindrical box with the vanishing 
boundary condition. The HFB equation is solved with 
the cutoff at the quasiparticle energy of 20 MeV in order 
to avoid huge test calculations in terms of computational 
amount. The canonical quasiparticle basis [17| is used for 
the general quasiparticles in the formulation mentioned 
in Sec. [Til an d the basis wave functions are constructed 
so as to include the unbound components according to 
the method of Ref. (67} after the HFB equation is solved. 
The unbound components are important for accurately 
obtaining the wave functions of the lesser occupied levels. 

Subsequently, trimming of the basis space is carried 
out by removing a small number of the canonical quasi- 
particle states with the least occupation probabilities in 
each (q, 7r, j z )-subspace so that the dimension of the sub- 
space is the same between the bases of the two nuclei if 
the corresponding original dimensions arc different. This 
process is necessary in our calculation, because the di- 
mension of each subspace is not a direct input to the 
HFB calculation, but the dimension is controlled by the 
cutoff quasiparticle energy. 

After this adjustment, the matrix elements of the uni- 
tary transformation are calculated according to Eqs. ([32]) 
and (|33p . The two canonical quasiparticle spaces are not 
identical in the coordinate calculation if the correspond- 
ing dimensions are identical, because the truncated space 
is determined sclf-consistently by solving the HFB equa- 
tion. Thus, a small correction is made to T IF1 and T IF2 
in such a way that the canonical quasiparticle wave func- 
tions associated with the state / obtained by the trans- 
formation (|31[) arc orthonormalizcd. The states with less 
occupation are mainly modified in this orthonormaliza- 
tion. The canonical quasiparticle wave functions are used 
for calculating T IF1 and T IF2 and the interaction matrix 
elements in the QRPA equation. 

Subsequently, we calculate the matrix D [Eq. (|35[) ] and 
the normalization factor of the HFB state [Eq. ([28]) ]. The 
D matrix is singular, when the two HFB states are or- 
thogonal as is seen in Eq. (|28[) . that is, the two nuclei 
share the same single-particle basis but have different 
configurations. This does not occur, however, as long 
as the HFB equations of the two nuclei are solved self- 
consistcntly. 



We solve the QRPA equation in the so-called ma- 
trix formulation [l7l . [67| . The two-canonical-quasiparticle 
spaces used in the QRPA calculation are not truncated 
in the test calculations after the trimming. This treat- 
ment enables clear separation of the spurious states asso- 
ciated with the particle number. We reduce the size of the 
canonical-quasiparticlc spaces using the tight cutoff of 20 
MeV so that all allowed combinations of the two canon- 
ical quasiparticle states are easily tractable. Thus, the 
discussion of the nuclear properties is out of the scope of 
this paper. Using the QRPA solutions, one can calculate 
the matrix elements of the generator of the QRPA ground 
state [Eq. ((26)) ] and the associated normalization factors 
[Eqs. (HOP and gTJ)]. The matrix X^'^ 1 or X^ K '^ F 
would be singular, if all of the forward amplitudes of a 
QRPA solution vanish, or a two-canonical-quasiparticlc 
component is not used in any of the QRPA solutions. 
However, this does not physically occur. 

The next step is the calculation of the contractions 
d62j) (J64J) using Eqs. ([32]), ([33]), and ([35]). The general- 
ized expectation values of high order with respect to the 
fermion operators arc calculated by generating the list 
of the indices of the canonical quasiparticle states used 
in the contractions systematically and recursively from 
a low order (refer to the proof of Wick's theorem jo8j). 
Finally, the overlap matrix elements ([3"8]) are calculated 
by using Eqs. ([ST]) - ([60]) . 



Truncation scheme 



Feature of parallel computation affects the answer of 
a question of what approximation is efficient. Equa- 
tion ([ST]) can be calculated by multiplication of three 
matrices having matrix elements X F " 1 * , Xjfy and 

(f\a F a^a^aj^i). These matrices are partitioned, dis- 
tributed to the cores of the computer, and handled by 
ScaLAPACK [69[ in the process of the multiplication. 
However, this approach is not efficient in the calculation 
of Eqs. ([52]) — ([60]) . This is because the redistribution of 
the matrix elements of C ( - KlV ^ F * and C^ 1 ^ 1 ^ 1 is neces- 
sary between the cores before the matrix multiplication is 
carried out using ScaLAPACK when different terms are 
calculated. This requires large computation times if the 
data size is large. Thus, we calculate the un-normalized 
overlap of Eqs. ([52]) — ([60]) truncating the two-canonical- 
quasiparticle states used in the summations without us- 
ing ScaLAPACK (parallel computation is still used). The 
efficiency of this truncation is high due to the reason for 
the u-expansion discussed in Sec.[lTj that is, \I) and \F) 
have different configurations at the Fermi surface. This 
efficiency is confirmed numerically in Sec. [V] Obviously, 
this approximation also holds good for Eq. ([ST]) . How- 
ever, that term is calculated without this approximation 
because the matrix multiplication using ScaLAPACK is 
very efficient. We introduce another independent trunca- 
tion of the two-canonical-quasiparticle states for the cal- 
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dilation of ([55)) because those states {(/xz/)} that 

are sufficient for the lower-order terms (|52|) and (|53|) are 
too many to handle for the six-fold summation with re- 
spect to these states in the calculation of F^ m , . 

Two truncations are used regarding Kir; one is that 
in the summations of Eqs. (|52j) — ([5(7)1 . and another is 
that of A4 [Eq. ([39]) ]. These truncations are treated in- 
dependently. The last one is the truncation with respect 
to v { f l7Tl)f and 4 Kl7ri \ as is shown in Eq. ((35J). Wc 
calculate the normalization factor up to an order higher 
than that of the un-normalized overlap, as mentioned in 
Sec. [HI Thus, we use six truncations in the calculation of 
the overlap matrix after the canonical-quasiparticlc bases 
are determined. These truncations are investigated nu- 
merically in Sec. [V] 



It is our intention to test two fairly different cases; in 
one case, the two nuclei have similar properties except 
for the difference in the proton and neutron numbers, 
and in another case, the two nuclei have fairly different 
properties. This difference can siginificantly affect the 
overlap matrix, because the matrix is not unitary, i.e., 
there is no normalization of the matrix elements. 

It is a physical feature of the region around 26 Mg that 
the sign of the ground-state quadrupole deformation is 
sensitive to the input parameters. We confirmed that 
the HFB ground states were axially symmetric using a 
three-dimensional HF-plus-BCS code as long as SkM* is 
used. 



C. Properties of test states 

We discuss the physical properties of the test states 
that we use in this study. 26 Mg and 26 Si are used for \i) 
and /) (If 1 )), respectively, with the Skyrme param- 
eter set SkM* [701 an d the volume pairing density func- 
tional [7l[ . Two sets of test states are used with different 
pairing strengths. The properties of the HFB ground 
states are shown in Table [1] and the pairing strengths 
are given in Table [ill The total dimension of the HFB 
space is ~330 including those with negative j z . 



TABLE I: Properties of HFB ground states of 26 Mg and 26 Si 
for test sets I and II. /3 P and A p denote the quadrupole defor- 
mation and the averaged pairing gap of the protons. Those 
with the suffix n correspond to the quadrupole deformation 
and the averaged pairing gap of the neutrons. 



Nucleus 


Pp 


A p (MeV) 


A. 


A n (MeV) 


Test set I 


26 Mg 


-0.199 


0.794 


-0.195 


1.510 


26 Si 


-0.224 


0.865 


-0.206 


1.402 


Test set II 


26 Mg 


-0.228 


0.779 


-0.234 


<0.001 


26 Si 


0.251 


0.011 


0.316 


0.259 


TABLE II: Strengths of pairing enei 
used. A cutoff quasiparticle energy 
HFB calculations. 


gy functional G p and G n 
of 20 MeV is used in the 


Nucleus 




(MeV fm 3 ) 




(MeV fm 3 ) 


Test set I 


26 Mg 




-150.0 




-270.0 


26 Si 




-270.0 




-200.0 


Test set II 


26 Mg 




-150.0 




-150.0 


26 Si 
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FIG. 2: (Color online) Ten diagonal overlap matrix elements 
having largest absolute values as functions of %If + 91/. 
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FIG. 3: (Color online) Ten diagonal overlap matrix ele- 
ments having largest absolute values as functions of ((Ol^ 1 + 

?nF 1 )/2) 6 . 
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|K| 

FIG. 4: (Color online) Convergence of 1/M, Eq. ([39]), with 
respect to \K\. Each value includes the contributions of both 
tv = + and — , and the terms with K < are also included. 
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(Color online) Energies of K* = + QRPA excited 
26 Mg and 26 Si in order corresponding to data in 
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FIG. 5: (Color online) Twenty largest absolute values of 
diagonal overlap matrix elements. Those up to the second 
order with respect to Cj 0+ ' and v^ + ^ are shown in descending 
order. The terms with (Knn) / (0+) of Eq. (J3H1) are not 
included. We used 91 F + 91/ = 134 and ((9T? 1 + f ttf 1 )/2) (i = 
3 x 10 s (see Figs. and 0. 



V. NUMERICAL TEST OF TRUNCATIONS 
A. Test set I 

We examine the effects of the various truncations sep- 
arately using test set I with (Ktt) = (0+). The conver- 
gence with respect to the number of the two-canonical- 
quasiparticle states used in Eqs. (|52j) and ((53]) is ex- 
amined without the second-order term with respect to 
v\ K ^ and v^p^ [Eq. ([54]) ]. The truncation of the two- 
canonical-quasiparticle states for F^ m , [Eq. ([55]) ] is in- 
vestigated by suppressing {K\tt{) ^ (0+) in the un- 
normalized overlap [-F^ m ,, Eq. ([60]) . is omitted for sim- 
plicity] . The effect of Eq. ([51)1 is also investigated using 
only (Kitti) = (0+) in the un-normalized overlap. On 
the other hand, when the terms with (Kitti) ^ (0+) arc 



FIG. 7: (Color online) Summation of squared backward am- 
plitudes of QRPA excited states corresponding to Fig. [5] The 
definition of the symbols is the same as that of Fig. [6] 

included, Eq. ([M]) is omitted. 

Let 9Tf and 91/ be the number of two-canonical- 
quasiparticle states associated with \F) and \I) truncated 
for calculating Eqs. ([52]), (53J), and ([55 ]) -([50 ]) . Since 
Fmm'i Eq- ([55|l . has six- fold summations with respect to 
(fxv), it is truncated separately as mentioned before. We 
show the convergence of the overlap matrix elements with 
respect to Dt/r + 91/ in Fig. [21 91f and DTj are controlled 
in the numerical calculation by using a cutoff occupation 
probability of the canonical quasiparticle states, and the 
states with larger occupation probabilities than the cut- 
off are used. The occupation probability is defined by the 
norm of the lower component of the quasiparticle wave 
function, and in our calculation, it is equal to the occupa- 
tion probability of the canonical state. The same value 
of the cutoff is applied for the two bases, and we have 
9Tf ~ 91/. It is seen that 9V + 91/ = 350 is sufficient for 
the convergence. The total number without the trunca- 
tion is ~3300, and thus, this truncation is fairly efficient, 
as has been discussed before. 
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FIG. 8: (Color online) Twenty largest absolute values of 
diagonal matrix elements of overlap in descending order with 
all values of (|isTi|7Ti), that is (0+) — (4+), and those with only 
(0+) (not necessarily in descending order), see Eq. (|38|l . The 
converged M value is used for both calculations. Negative 
parity contributions are not included (see text). The QRPA 
state m = 1 (2) corresponds to m = 2 (1) of Fig. [5] 
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FIG. 9: (Color online) Twenty largest absolute values of 
diagonal matrix elements of overlap in descending order for 
QRPA states of (Kn) = (2+) with only zeroth-order terms 

fn i \ (2 + ) 

with respect to v\ and v F and those also including first- 
order terms. The terms with (Ki'Ki) ^ (2+) in Eq. (|38p are 
not included. Truncation was made at VIf + 91/ = 350. 



Figure [3] illustrates the convergence of the diagonal 
overlap matrix elements with respect to the number of 
the two-canonical-quasiparticle states used for the cal- 
culation of i 7 ' T J lm , (|55|) . The term DTj^ 1 is the number 
of those states associated with \F), and Vtf 1 is that as- 
sociated with \I). The terms ((D^ 1 + 9tf 1 )/2) 6 is the 
number of the terms of the six-fold summation with re- 
spect to (fiv) of Eq. dHJ [only (Kitti) = (0+)]. A 
fairly stable convergence is obtained, and the value of 
((OTp +mf 1 )/2) 6 = 5 x 10 7 is sufficient for convergence. 
This implies that OT^ 1 is at most 20 with D^ 1 ~ Otf 1 , 
and thus, the second-order term (|54p can be considered 
as negligible. 



FIG. 10: (Color online) The same as Kn = 0+ term in Fig. [8] 
but for the terms proportional to q {f\i) q , q denoting proton 
or neutron, not included. 



The normalization factor 1 j M. = NiMf does not have 
a mechanism of fast convergence with respect to 91f +91/ 
unlike the un- normalized overlap matrix element, as men- 
tioned before. Indeed, we found that no truncation was 
possible to satisfy Eq. (|7ip. thus, 1/A4 was calculated 
without that truncation. Figure [4] depicts the conver- 
gence of 1/M with respect to \K\, thereby indicating 
that \K\ up to 3 is sufficient. Further, we found that the- 
normalization term up to the second order in Eq. (1401) . 



1 



(72) 



Kim 



was 3.843, and that for |/) was 4.053. The fourth-order 
term for \i) [Eq. (|l8"]) ] was 0.980, and that for |/) was 
0.838. The first term of Eqs. (gHJ) and (@9j) (called as the 
unlinked term) was found to be larger than the second 
term by a factor 2—3. Thus, order estimation is possible 
for the normalization term of the sixth order 



36 ^ 

K 171-1 



»| (fl^M) 3 -)) 3 | i>) (73) 



which is not included in the calculations of the overlap, 
by considering the unlinked terms included in Eq. (|73p 



i y: (Tr ( c(— )c^-)t)) 3 x ( 4 6 8 



6, K^0, 
K = 0. 



(74) 



Equation d74j) gives a value of 0.088 for \i) and 0.113 for 
|/). We ignore this order of contributions. The exchange 
terms in the second-order normalization terms (|44D and 
(|47)) were found to be around -0.13 for (Kn) = (0+), 
and all of the absolute values of the terms with other 
(Ktt) were smaller than 0.01 with the tendency that the 
larger was the value of \K\, the smaller were the absolute 
values. These terms are also negligible. 

The major diagonal overlap matrix elements are shown 
in Fig. It is seen that the contribution of the second- 
order term (|54[) is negligible, and that of the first-order 
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terms ([52]) and ((53)) arc not significant for the small ma- 
trix elements. The zeroth-order term (|51|) is sufficient in 
most of the matrix elements omitted in that figure. 

Figure [S] shows the energies of the QRPA excited states 
with (Kit) = (0+) of 26 Mg and 26 Si. It is seen from 
Figs. [5] and O that the major diagonal matrix elements 
of the overlap arise from the states with energy lower 
than 15 MeV. The charge symmetry of the two nuclei 
is also obvious — it is perfect in the low-energy region. 
Figure [7] shows the summations of the squared backward 
amplitudes of the QRPA solutions. These energies and 
backward amplitudes indicate that the states m = 1 and 
2 are the spurious solutions associated with the particle 
number. 

TABLE III: Contribution of (|Ki|7Ti) to first-order term [ 
Eq. (J52J) plus Eq. l[g3]l ] for m = 1. see Fig. El The sum- 
mation of these contributions is equal to 0.32 x 10 -2 , and the 
summation except for (Kvk\) = (0+) is 6.24 x 10 -2 . 



\Kx\-kx 


Contribution to the first-order term 
(xl(T 2 ) 


0+ 


-5.93 


1+ 


3.26 


2+ 


0.47 


3+ 


0.63 


4+ 


1.90 


5+ 


-0.002 



The contribution of (Kitti) ^ (0+) to the first-order 
terms ([52]) and ([53]) of the major overlap matrix elements 
are shown in Fig. [51 We calculated that of (Kitti) = (0— ) 
and (1 — ) and found that it was smaller than that of 
the positive parity by at least an order of magnitude, 

I 



1 1\ F F I\ It I -\ 



The -ftf-quantum number of \f) q and \i) q (q = p or n) 
is zero, thus, the first term of Eq. ([77]) vanishes for K 
values of the QRPA state other than zero. Therefore, the 
overlap matrix elements of Fig. [5] arise from the break in 
the particle-number conservation. 

Subsequently, we calculated the overlap matrix of the 
(Kir) = (0+) states without including the terms propor- 
tional to q{f\i)q, and the corresponding results arc shown 
in Fig. 1101 Upon comparing the largest overlaps except 
for those of the spurious states observed in Fig. [TO] and 



thus, only the positive parity is used. The contribution 
of (Kitti) ^ (0+) is very small to all the diagonal matrix 
elements except for those of the spurious states. This 
extreme sensitivity of the overlap of the spurious states 
provides us with one more reason why our method should 
be applied only to the cases for which the break in the 
particle number conservation is small. Table IIIII sum- 
marizes the details of the contribution to the most sensi- 
tive matrix element. The (|i£f |7r)-dependence is irregular; 
however, eventually the contribution becomes negligible 
for the large value of \K\. We also examined the contri- 
butions of the quasi-boson terms and the exchange terms 
of {flv^^OmOUf). The absolute value of the quasi- 
boson term is larger than that of the exchange term by a 
factor 2—10 in many of the major matrix elements. Thus, 
the quasi-boson term is the leading term. 

We also calculated the overlap matrix elements of the 
(Kit) = (2+) states (Fig. [SJ). The result of the zeroth- 
plus-first order term can be compared with that of the 
(Kir) = (0+) term of Fig. [8j The several largest values 
of (Ktt) = (2+) are 50—60 % larger than those of the 
real states of (Ktt) = (0+), and the small ones in the 
tail of the curve are comparable. We discuss the origin 
of the overlap of the (Ktt) ^ (0+) states based on the 
zeroth-order term (|5Tj) . The HFB wave function can be 
expressed as a direct product of the proton and neutron 
wave functions, that is, 

l/> = l/> P ®l/>«, (75) 
|») = \i) P ® \i)n, (76) 

where \f) p (\i) p ) and /}„ (\i) n ) denote the proton and 
neutron wave functions. Since \x and v (// and v') in 
Eq. ([ST]) are like particles, the generalized expectation 
value used in that equation is written as 



fiv: protons and 
fx'i/': neutrons, 

\iv\j! v'\ protons, (77) 



I 

those corresponding to the result labeled (Ktt) = (0+) 
in Fig. [8j it is observed that 55% of the overlaps of the 
(Kn) = (0+) states arise from the equations that do 
not vanish in the case where the particle number is con- 
served. However, it is to be noted that removing the 
terms proportional to q (f\i) q artificially is not justified 
because the consistency between the equations derived is 
ignored. The optimal approach is to carry out the parti- 
cle number projection of the many-body wave functions; 
however, this is out of the scope of this paper. 



' p(f\a?a%\i) P n(f\aj,aj\i) n , 

pifla^a^ajajl^p n {f\i) n , 
the same equations but with 
, the protons and neutrons exchanged. 
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FIG. 11: (Color online) Ten diagonal overlap matrix ele- 
ments having largest absolute values as functions of 91 f + 9Ti 
obtained using test set II. The terms with (K\ir\) 7^ (0+) 
are included in the first-order terms (|52[) and (|53p . and the 
second-order term (|54|l is not included. 
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FIG. 12: (Color online) Twenty largest absolute values of 
diagonal overlap matrix elements obtained using test set II. 
The terms with (-Ki7ri) 7^ (0+) are not included in the first- [ 
Eqs. (|52[) and (153[) ] and second-term (|54[) . and the condition 
01f + 01/ = 200 is used. 



B. Test set II 

In this subsection, we discuss the results obtained using 
test set II (Table [J). Figure [Til illustrates the m F + 91/ - 
dependence of the major overlap matrix elements. The 
convergence is slow compared to that of test set I shown 
in Fig. [21 The possible reason is that the canonical quasi- 
particlc basis is rather different between the two nuclei, as 
is seen in Table [H so that the basis wave functions with 
low occupation probability have certain components of 
the wave functions with larger occupation probabilities 
in another basis. 

Figure [T^] shows the major diagonal matrix elements 
of the overlap. A couple of siginificant differences from 
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FIG. 13: (Color online) Twenty largest absolute values of 
diagonal overlap matrix elements obtained with and without 
(A"itti) = (l+)-(5+) using test set II. The condition 91f + 91j 
— 200 is used, and the second-order term (1541) is not included. 
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FIG. 14: (Color online) Energies of K* = 0+ QRPA excited 
states of test set II. For the definition of the symbols, see 
Fig. El 



Fig. [5] can be observed; one is that the values in Fig. [12] 
are an order of magnitude smaller than those in Fig. [5] 
This implies that if the two ground states are quite dif- 
ferent, the overlap matrix elements become very small. 
The values of S defined by Eq. (J66J) arc ~1CT 2 (test set I) 
and ~1CT 5 (test set II) with dim(G) ~ 1650. The second- 
order term ([54)) and the terms with {K\-k\) 7^ (0+) in the 
first-order terms (f52jl and ([53]) are not included in these 
calculations. The value of S confirms that the all the 
matrix elements of the overlap are reduced significantly 
when the difference in the structure of the ground states 
of the two nuclei increases. 

It is noteworthy to compare the above values with that 
of the simple model discussed in Sec. IIIII The S value of 
that model is 2/1650 ~ 10~ 3 , and this value is located 
between those of test sets I and II. This order implies that 
the ground states of the two nuclei of test set I are closer 
to each other than those of the simple model, and those 
of test set II differ considerably from those of the simple 
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FIG. 15: (Color online) Summation of squared backward 
amplitudes of QRPA excited states corresponding to Fig. 1141 
For the definition of the symbols, see Fig. [6] 

model. The probable reason for the former relation is 
the singlc-particlc-configuration mixing arising from the 
similar pair fields of test set I (some components of the 
nuclear wave functions should be shared), and that for 
the latter relation is the quite different mean fields of test 
set II (the order of the single-particle levels is different). 
A crude guideline is given for the overlap matrix from 
this argument as 

- log(dimG) - 2 < log 5 < - log(dimG) + 1. (78) 

Another difference from the result of test set I is that 
the zeroth-order term (|5T|) is sufficient in many of the 
overlap matrix elements. The energies and the summa- 
tion of the squared backward amplitudes of the QRPA 
solutions are shown in Figs.ll4land[T5l respectively. The 
results indicate that the QRPA solutions m = 1 and 11 
are the spurious states associated with the particle num- 
ber. Upon comparing Figs. 171 and[l"5l we observe that the 
summation of the squared backward amplitudes of test 
set II is smaller on average than that of test set I except 
for the spurious states. Thus, the difference in the effect 
of v^ 1 ^ 1 ^ and v^ Kl7Tl ^ can be explained by the backward 
amplitudes [see Eq. ([26]) ]. 

VI. SUMMARY 

The overlap matrix elements of the QRPA states based 
on the ground states of different nuclei have been calcu- 
lated using the QRPA ground states explicitly. Our idea 
for handling the QRPA ground state is to expand it with 
respect to the generator — this approach appears to be 
a highly feasible method, and the expansion is proba- 
bly the only feasible method. Further, certain analyti- 
cal properties of the overlap matrix have been discussed. 
The non-unitarity is the exotic mathematical property of 
the overlap matrix in the QRPA. The truncation scheme 
used is explained in relation to the parallel computation, 



and the calculations are performed using relatively light 
nuclei with the Skyrme and the contact volume pairing 
energy functionals. The truncations of the calculation 
have been examined carefully and justified numerically. 

The computation provided the three following bene- 
fits: 1) The truncation of the two-canonical-quasiparticle 
space is efficient in the calculation of the un-normalizcd 
overlap matrix. The normalization factor requires cal- 
culation with no truncation of the two-canonical quasi- 
particle space used in the QRPA calculation, and on the 
other hand, the calculation is reduced tremendously by 
using the identical bra and ket states, 2) the inclusion 
of up to the linear term with respect to the generators 
of the QRPA ground state is sufficient in the expan- 
sion of the un-normalizcd matrix elements, and 3) the 
(A'i7Ti) ^ (Kit) terms contribute negligibly to most of 
the un-normalized overlap matrix elements. The reason 
for the first benefit is obviously independent of nuclei, 
and the second benefit should be applicable to any nuclei 
for which the QRPA is a good approximation. 

As for the normalization factor, the maximum value of 
\K\ = 3 is sufficient in our test calculations. The terms 
up to the fourth-order with respect to {jj ifl7ri - ) and v^ l7r '^ 
were calculated in Aff and J\fp, and the next-order terms 
were estimated to be negligible. 

Certain selection rules on the quantum numbers of the 
two-canonical quasiparticle states are used in the various 
terms, and the terms with (A'itti) ^ (Kir) are subject to 
one more condition than the terms with (A'itti) = (Kir). 
Because of this difference, the former case has less terms 
than the latter case. Thus, the terms of the former case 
seem less coherent, and this is the only explanation as 
regards the third benefit. Therefore, if the above expla- 
nation is correct, all of these benefits should also hold 
for heavier nuclei. Considering these advantages and the 
recent development of powerful parallel computers, there 
is no reason to avoid performing calculations using the 
explicit QRPA ground state if accurate calculation is nec- 
essary. Two sets of nuclear wave functions were used, 
and it has been shown that the overlap was sensitive to 
the difference in the wave functions of the initial and fi- 
nal states. The feasibility of the calculation has been 
demonstrated in both the test cases. 

We have included as many terms as possible in our 
calculations. As a result of this manner, certain non- 
quasi-boson (exchange) terms are included, while others 
are not. The generators of the QRPA ground state were 
obtained by the quasi-boson approximation. It does not 
seem possible to extract the matrix element Cj^™?^ 1 
in the isolated form unless the exchange terms are ig- 
nored. According to our experience of the calculation 
of the equations including both the quasi-boson and ex- 
change terms, the exchange terms are not as significant 
as the quasi-boson terms perhaps because of a reason 
similar to that for the third benefit discussed above. 

The code is developed in such a manner that memory- 
shortage issues do not occur if applied to heavy nuclei, 
and the parallelization efficiency is good. Thus, the ap- 
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plicability of our method to heavy nuclei is a matter 
of availability of core-hours. We are preparing to ap- 
ply our method to the calculation of the nuclear matrix 
elements including the phase-space factor for a dozen of 
the 0v(3(3 decays. Finally, it should be possible to apply 
our method to calculate the overlap to the pn-QRPA. 
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